Dipolar Self Assembly in Coherently Fluctuating Fields 



Igor M. Kulic^ and Miodrag L. Kulic^ 
^^CNRS, Institute Charles Sadron, 23 rue du Loess BP 84047, 67034 Strasbourg, France and 

Goethe-University D-60438 Frankfurt am Main, Germany 
(Dated: April 23, 2013) 



Institute for Theoretical Physics, 



o 

(N 

u 

Oh! 

< 

(N 



o 



-a 

c 
o 
o 



> 

o 
in 

o 



x 

S3 



When two polarizable objects interact in spatially randomly fluctuating fields, we speak of Van 
der Waals-like forces. But what happens if the random excitation becomes spatially coherent? 
We show that such a spatially coherent fluctuation interaction (scFI) displays remarkable physical 
features distinguishing it from its sister, the incoherent fluctuation interaction (icFI). In particular, 
we investigate how scFI gives rise to complex, hierarchically structured, self-assembled materials, in 
contrast to the usual phase-separated droplets generated by icFI. We discuss a recent example of 
scFI, implemented with superparamagnetic colloids in triaxial magnetic fields by Martin et al. and 
explore the main steps in the hierarchical formation process of coherent colloid self assemblies, from 
chains, through membranes, to large-scale foamy superstructures of the material. We predict that 
such coherent "dipolar doughs" exert positive pressure, can easily swell against gravity on cm-scales 
and we derive their equation of state. 



PACS numbers: 82.70. Dd, 81.16.Dn, 82.70.Rr 

The basic motif behind the formation of most solids 
and liquids is the interplay of a long-range attractive 
and a short-range repulsive force. As happens often in 
physical systems, simply reversing the signs of the in- 
teractions is a real game changer. A competing short- 
range attractive and long-range repulsive force can give 
rise to profoundly different, complex, sponge-like struc- 
tures ranging from gyroid phases in block copolymers 
[l|, labyrinthine phases in ferrofiuids Q to nuclear pasta 
phases in gravity-collapsed neutron stars [3| . By impos- 
ing spatially coherent conically rotating (triaxial) mag- 
netic fields Martin and coworkers have discovered [J, [5[ 
that superparamagnetic Ni- colloids, placed in a finite- 
volume container, self-assemble into unusual foam-like 
structures, displaying a local cohesive force on small 
scales but with an overall large-scale tendency to re- 
pel and rise against gravity, much like a swelling dough 
(cf. Fig. 4a below). Owing to the complexity of the 
spatially coherent oscillating field-induced magnetic in- 
teraction and its inherent anisotropic dipolar nature, 
the structures formed on the smallest scales unexpect- 
edly resemble mono-layered colloidal membranes [4j-[6| 
- in sharp contrast to the chains predominating in the 
usual time-constant fields. After a continued growth 
phase the colloidal membranes begin to merge into larger- 
scale foam-like superstructures spanning the whole sam- 
ple container. What Martin et al. have discovered in 
their seminal work 4, 5], is a particular instance of a spa- 
tially coherent fluctuation interaction (scFI): a spatially 
homogeneous dipole-exciting field with certain statistical 
properties (its square is statistically isotropic) varies in 
time and gives rise to an effective interaction (cf. Fig.l) 
of amazing complexity with a significant potential for 
micro- and nano-scale self-assembly. 

This letter aims at bringing order into our theoretical 
understanding of scFI by outlining a certain simplicity 
underlying its seeming complexity. We will explore how 



scFI's intrinsically strong many-body interactions drive 
the formation of chains and membranes and finally how 
complex membrane-membrane pair interactions give rise 
to swelling " dipolar doughs" , whose equation of state we 
derive at the end of the letter. 

1. Coherent vs. Incoherent Fluctuation Interaction. 
How can the scFI (cf. Fig.l) generate such complex 
structures while its sister - the incoherent fluctuation in- 
teraction (icFI)-merely forms phase-separated lumps or 
droplets of matter within a two phase system [7|? To 
answer this question we consider dipolar magnetic (di- 
electric) particles without permanent moments [8|], such 
as colloids, which are placed in a temporally fluctuating 
magnetic (electric) field Bo,; = jUoHo,i(t), which changes 
on a typical timescale t u . The latter fulfills the condi- 
tion tm < t u < T n with tm the dipolar relaxation time 
and t v the characteristic time for the colloid's motion 
in the surrounding viscous fluid. Under these conditions 
the beads' magnetization is equilibrated, while their cen- 
ters of mass respond much more slowly and feel a net 
time-averaged force. The free energy functional T for in- 
teracting beads with magnetization M; in imposed fields 
Ho,j is then given by [9(] 

—T1&LH H., ) = Y. (^ M «M lM - - M*Ho,i) Vi (1) 



\ J2 MAjMjVi, 



where the summation goes over N beads, fig is the vac- 
uum permeability, Vi is the volume of the i-th bead 
and Xb is the bead susceptibility tensor. For an el- 
lipsoidal colloid (or colloid aggregate), one has \b — 
(1 + LbXmb) Xmbt where Xmb is the colloid's material 
susceptibility - which can take very large values in multi- 
domain (or superparamagnetic) beads - and Lb is its de- 
magnetization tensor. The first term in Eq|T] comes from 



the ferromagnetic domain entropy within beads, while 
the bead-bead magnetic interaction is contained in the 
third term [l(|. Here we keep only the magnetic dipole- 
dipole interaction [I, Q, given by Ty — y>y£(by), with 
tfij = (p(Hi, Hj) — Vj/Air |Rjj | , Ry = Rj — Hj and the 
tensor i(by) = 1 — 3by (8) by (dyadic product) with the 
"bonding" unit vector by = Ry/ |Ry|- After minimiz- 
ing T with respect to Mj for identical beads {V\ = V&), 
followed by time-averaging over a spatially coherent ex- 
citation, H ,i = H , H^QaHjflf} = 5 a pH$ (i,j = 1, .., JV; 
a,/3 = x,y,z), one obtains a physically transparent ex- 
pression for the effective energy in the scFI case 
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FIG. 1: a) The standard incoherent (icFI) and the b) spa- 
tially coherent fluctuation interaction (scFI) are both induced 
by field fluctuations - but with different spatial correlations. 
c)+d) While the 2-body forces appear superficially similar in 
both cases, the 3-body forces have different angular character 
and are longer ranged for scFI, Eq.©. 



spatial trace of the (non-local) effective susceptibility 

Xeffnj = Xb (Ri| {i+Xbf)- 1 |Rj> With (R,| f |Rj) = fy. 

As Xeff,ij describes the coupling of the i-th and j-th bead 
in presence of all other beads, it gives rise to many-body 
effects. These are crucial for the hierarchical structuring 
of colloids, as confirmed in numerical simulations [J, l5| 
and experiments [4j-|6|, where dimers, molecules, chains, 
branched chains, membranes and doughs are formed, in- 
stead of droplets and close-packed 3D crystal structures. 
The elegant "trace- formula" Eqj2] forms the basis for our 
further study of scFI systems. Note here that, in con- 
trast to scFI, in standard incoherently excited FI (icFI) 
systems, with -H^oa-H^o/a — S^Su^Hq (note the 6ij), the 
summation of Tr {x,eff,ij} hi T includes only i — j terms, 
giving rise to the usual Van der Waals(VdW) -like forces, 
preferring the formation of simple droplets only pj. In 
the following, we study first the bead-bead interaction, 
then the formation of chains and membranes, and finally 
the assembly of membranes in the form of 3D doughs. 

2. Dimer formation. In a first step, let us con- 
sider a very dilute system, where the pairwise bead- 
bead interactions should (naively) dominate in J- sc di- 
For two spherical beads (with indices 1,2), we have 

L a fi = (l/3)5 a /3, (Xb)a/3 = Xb0~a/3 (with Xb < 3) 

and Eq.©, gives Tr {jgj ' /12 } = 3 X6 (1 - X b^b)/(l + 

X&¥>&)(1 ~ 2xb^h), where tp b = V&/47r|Ri 2 | . This means 

that the bead-bead interaction is attractive for \b > 0. 

For large distances |R 12 | > V 6 1/3 , the interaction is e h2,3 cf. Fig.ld). Remarkably, F^l is of the same 

jffi = -(3/87r 2 )(VbXb) 3 i? 2 |Ri2r 6 giving rise to an order in R and Xb as the 2-body interaction P s l% 



experiments [4|-[6j , which show a clear tendency for chain 
and membrane formation. What is the microscopic origin 
of these complex structures in scFI systems? The reason 
hides in the specificity and the strength of the many- 
body interactions, which mark the real difference between 
scFI and icFI systems. In a typical icFI system, the 3- 
body interaction is described by an Axilrod- Teller type 3- 
body potential P^-r = C123/ IR12I 3 IR23I 3 |R.3i| 3 , with 



C123 oc xl (1 + cos ^i cos #2 cos #3) [ll|. Due to its weaker 
~ x 3 R 9 scaling, it is typically small and overridden 
by the 2-body VdW-like interaction ~ X 2 -R~ 6 [l2|, giv- 
ing rise to close packed droplets in icFI systems. In 

— (123) 
sharp contrast, the 3-body interaction Fg cFI in scFI sys- 
tems scales very differently with distance, as seen from a 



O 



xM) 



expansion of Eq. ([2]) |lj 



r. 



(123) 



scFI 



3 COS 2 gi - 1 

32tt 2 ^ u " 6 ^^ 



2„,3t/-3 



rl*oH XbVb 



|Ri2|°|R 



perm. (3) 



13 



where we sum over all 3-index permutations (for angles 
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isotropic and short-range attractive VdW-like attractive 
force F12 = —(dT{ 2 /<9Ri) <x IR12P ■ This interesting 
two-body result was first obtained in the seminal papers 
[J, l5| and confirmed experimentally [6] . 

3. Many-body effects. At first glance, the isotropic 
VdW-like two-body interaction oc 1/ IR12P appears to 
favor droplet-like, bulk structures. However, this contra- 
dicts numerical simulations by Martin et al. [4j, |5| and 



-Xb/ IR12I ■ In fact, the 2-body interaction is formally 
contained in Fg C pi (indices 1^2=3) showing us the pit- 
fall in the "Dimer formation" section: For scFI the 2- 
body interactions are physically inseparable from the 3- 
body ones - at any x&! Interestingly, ^ cFI has a sim- 
ple angular dependence explaining scFFs tendency to de- 
stroy 3D bulk structures: the — cos 2 9% term favors 0. L = 
or 7r - i.e. the colloidal chains found in the experiments 
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FIG. 2: Left panel: scFI leads to the formation of colloidal 
clusters with random orientations growing over time [J, Q. 
At the early stages, smaller clusters are born as linear chains 
and transformed into membranes beyond a critical size N c (x). 
Right panel: Phase diagram of clusters of size N and suscep- 
tibility X- 



4- Chains and membranes. From the 3-body forces of 
EqJ3]we understand now intuitively why chains form ini- 
tially. To capture their transition to membranes, higher 
O (xifi) terms beyond Eq[3]are necessary. It is however 
conceptually more instructive to take a more macroscopic 
approach, where chains/membranes are replaced by pro- 
late/oblate ellipsoids. Here, Xeffnj m Eq[5] is replaced 
by its shape-dependent continuum limit x , which is 
spatially homogeneous inside ellipsoidal samples with di- 
mensions a = b ^ c and volume V = (4-7r/3)a 2 c. The 
free-energy T S eif and x are in this case given by 
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where L is now the demagnetization tensor of the com- 
posite ellipsoidal structure, with TrL = L a + L + L c = 1 
and Li its eigenvalues [9(. Here, x ~ Xfc/(1 — X&//3) -1 
is the composite material susceptibility, which is clue to 
local field effects in the dense bead aggregate, with /(<1) 
the volume fraction of beads within the membrane/chain. 
In the continuum limit, when the number of beads 
is large (N ^> 1), chains can be considered as ex- 
treme prolate ellipsoids with a long half-axis c ~ ND, 
a short half-axis a ~ D and the demagnetization factors 
L a = L b > L c = L ch ~ (\nN)N~ 2 . In that case one 
has Tr {x ch } = x(l + L^xT 1 + 2*[1 + (1 - L ch ) X /2]- 1 
in Eq|U Membranes, on the other hand, can be con- 
sidered as extreme oblate ellipsoids with two identical 
half-axes a = b ~ ^/ND , c ~ D and L c 3> L a — 
L b = L m ~ N- 1 / 2 . In this case one has Tr{ X m } = 
2x(l + L m xY x + x[l + (1 - 2L m )x] _1 in Eqg] In both 
cases, for chains and membranes, T S elf is dominated by 
the smallest demagnetization factors L m / ch — > 0, in the 
limit N — > oo. However, since for membranes two demag- 
netization factors vanish, while for chains only one van- 



ishes, the energy of a membrane is always smaller than 
that of the chain for N — > oo, i.e. T m < T ch '. For finite 
N, the chain-membrane transition is reached on the line 
N = N C ( X ), where P n {N c ) = T ch (N c ). As seen from 
Fig. 2, N c grows with the composite susceptibility x- For 
X ~ 3 (for "Dynabeads" with Xb ~ 1 M)i we estimate 
N c w 10— 15 - close to the experimental value N c ss 10 for 
the chain- membrane transition [6(. The agreement sug- 
gests that the continuum approach quantitatively cap- 
tures the behavior for finite N. 

5. Interaction of membranes - Once they emerge, what 
is the fate of the membranes as they continue grow- 
ing? How do they interact and mutually order during 
growth? To answer these questions, we need the 2- 
membrane interaction for arbitrary membrane orienta- 
tions 1X1,2 and anisotropic susceptibilities (of thin oblates) 
Xi = x(l + LiX) 1 ■ Up to linear order in tp m , one 

expands Xeff.i ~ Xi (1 - <Pm t(bi 2 )Xa ) ( and Xe//,2 
for 1— >2). For identical, distant membranes (y ml w = 

1/3 

V m , |Ri2| ^ Vm ), the O (<p m ) interaction energy in 
Eq.© reads 
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with Xmin / max the minimal/maximal eigenvalue of the 
membrane susceptibility tensor x , 7 = Xmin/Xmax and 
geometric factors C\ — ri! -bi 2 , c 2 = n 2 -b 12 , c 3 = iii -n 2 . 
This 2-membrane interaction is angle dependent and re- 
pulsive in many configurations - see Fig. 3. Notably, for 
fixed |Ri 2 | ,-Fint becomes minimal for the orthogonally 
"twisted" membrane orientation with iii _L n 2 , iii _L bi 2 
and n 2 _L bi 2 (ci/2/3 = 0). The "twisted" membranes 
attract each other since T^ < (up to the point of mu- 
tual contact), as in the " coplanar" case, yet the " twisted" 
configuration has lower energy. This interesting result 
should affect the kinetics of membrane formation: If two 
distant membranes start growing within a large L they 
will rotate to a 90° position before touching. Therefore, 
in finite £s, some type of glassy state in their orientation 
may be kinetically favored. In other relevant configura- 
tions, such as the " top" , with two out of plane parallel 
membranes (ci/2/3 = 1) or the "generic" one (cf. Fig. 3), 



the interaction is repulsive with < J-\ 
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<T\ 
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6. Formation of doughs and equation of state. Simula- 
tions and experiments [J, [5j provide empirical evidence 
for a hollow dough-like superstructure in finite Ls (cf. Fig 
4a) . What is the physical mechanism behind such " dipo- 
lar dough" formation? We have seen above that large 
aggregates prefer to form membranes, and that these 
membranes mutually interact. Specifically, when two dis- 
tant membranes are stacked over each other, they repel 



a 



00 



§0 



• /? 



12 




"Coplanar" 





"Generic" Repulsion 




Attraction 



"Twisted" Ground-State 




FIG. 3: The scFI is unexpectedly complex in various ge- 
ometries, a) The 2-bead interaction is simple in compari- 
son, always attractive and oc —R~ 6 in leading order, b) The 
bead-membrane interaction is longer ranged and anisotropic 
in relative coordinates, changing scaling to oc ±_R -3 . The 
interaction flips sign from attractive to repulsive on a cone 
with opening angle C rit ~ 54.7° [13J (c). The 2-membrane 
interaction is complex, attractive or repulsive depending on 
orientation, cf. Eq[5] with a ground state in the "twisted" 
configuration. 



FIG. 4: (a) The large scale structure of a dipolar dough (from 
0,13) in experiment (left) and in simulation (right), (b) The 
theoretical 3D shelf-model for dough's structure. 



ergy (for X » 1) ?int « (S/8Tr)V b tot f vX 2 Po, with 
po = -Bg/2/io and the lattice interaction constant [13| 



s = E^i ( c i + 4 + \ c l ~ c i c 2C3 - f ) Iom/RhI ~ 10 
obtained numerically from EqjSJ 

The self energies of finite membranes depend on their 
demagnetization factor L m (see discussion above), which 
for the shelf structure is given by fy as L m sa oc m fv for 
fy < 1 with a m « 1/4. For \ > 1 and L_ m x < 1 (dilute 
limit) the self-energy is given then by J^seif ~ const + 
Q-&X 2 fvPoVb"* ■ The corresponding contributions to the 
pressure are p se i f s=a 0.5f v x 2 Po and p int ~ OAf v x 2 po 
respectively, thus giving p m f v x 2 Po- The compressibil- 
ity of the dough is n = l/2p. For fy w 5 • 10~ 2 and 
Bo = 10 mT, we obtain p ~ 10 Pa and a large « ~ 0.1 

volume fraction in the container fy — V^ ot /V <C 1, with Par 1 . Since k oc i?^ 2 the compressibility is sensitive to 

V b * ot the total volume of all beads and V the container's 



{F^nt > 0)- I n the opposite limit (contact distance), a 
simple estimate implies their preference to split as well 
[14j . It is this remarkable reluctance of membranes to 
mutually stack that in fact sets the microscopic struc- 
ture of the " dough" . The dough is formed out of the 
thinnest possible membrane patches, whose thickness is 
collapsed onto the smallest available physical scale - the 
colloid size D. The characteristic lateral size om of these 
membrane patches ,on the other hand, is set by the bead 



volume. By assuming a cubic shelf structure as an ansatz 
(cf. Fig 4b) one obtains a patch size om = ZD/ fy. 

Since the membranes tend to grow, yet repel on the 
average, the container walls will feel a net pressure 
P (= -(dF/dV)) > [H, where T = T selj + T mt 
has contributions from both the membrane self-energy 
Fseif (from EqH]) and their pairwise interactions Tmt 
(Eq.(|3Jl). Since both attractive and repulsive mem- 
brane interactions are present, it is a difficult endeavor 
to find the exact equilibrium structure, even in nu- 
merical calculations. However, using the simple shelf- 
like cubic lattice ansatz with the mesh-size om (Flg-4), 
we can qualitatively mimic the dough structure [J, Q. 
For this regular lattice, we can sum the pairwise mem- 
brane interactions and obtain the total interaction en- 



the strength of excitation Bq. Finally, in a gravitational 
field the maximal height /i ma x of the dough is reached for 
p ~ Pgr, with the gravitational pressure p gr = Apghfy- 
For water immersed JVi-beads [J, |5( one has the density 
contrast Ap ~ 8 • 10 3 kg/m 3 and x ~ 10 m a densely 
packed membrane. For fy ~ 5 • 10~ 2 , B = 20 mT and 
X = 10 we predict that the dough rises significantly, with 
h° ~ 1 cm 0. 

In conclusion, we have studied, in the framework of the 
effective nonlocal susceptibility Xij j the formation of hier- 
archical superstructures in dipolar colloids driven by the 
spatially coherent fluctuation interaction (scFI). The pro- 
nounced many-body interactions (contained in Xij) gi ye 
rise to a growth of anisotropic assemblies - chains, then 
membranes once a critical cluster size N c (x) is reached. 
In a spherical container of finite size, smaller membrane 



patches of size ~ D/fv are formed, which, contrary to 
the case of attracting beads, repel on average, thus giv- 
ing rise to "dipolar dough" structures. The dough ex- 
erts a positive pressure on the walls of the container due 
to the tendency of membranes to increase their surface 
areas as well as their mutual repulsion. The "dipolar 
dough" represents an intriguing state of colloidal matter, 
formed by a subtle interplay of an attractive local interac- 
tion and a net repulsive longer range force. Remarkably, 
both types of forces are born out of a single, conceptu- 
ally simple interaction - the scFI (EqJ5]). The equation 
of state p ~ fyX 2 Bo gi ves r i se to a large compressibil- 
ity, highly sensitive to changes of magnetic field, open- 
ing doors for applications. Dielectric scFI analogues, like 
light-illuminated plasmonic nanoparticle systems, are ob- 
vious next challenges. Here, dynamic susceptibility ef- 
fects near resonance frequencies should add new physics 
needed to understand the shape-resonant self-assemblies 
found in Ref. [17| . The complexity of the quasi-static scFI 
considered here makes us wonder about the laws of the 
nano-cosmos emerging once dynamic effects are switched 
on. 
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